to appear as Rapid Com. in PRB (2010) 



o 

(N 
X> 

PL, 

in 

(N 



CO 



I 

O 

o 



> 

in 
o 
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High field electron spin resonance, nuclear magnetic resonance and magnetization studies ad- 
dressing the ground state of the quasi two-dimensional spin-1/2 honeycomb lattice compound 
InCu 2 /3Vi/303 are reported. Uncorrelated finite size structural domains occurring in the hon- 
eycomb planes are expected to inhibit long range magnetic order. Surprisingly, ESR data reveal 
the development of two collinear antiferromagnetic (AFM) sublattices below ~ 20 K whereas NMR 
results show the presence of the staggered internal field. Magnetization data evidence a spin reorien- 
tation transition at ~ 5.7 T. Quantum Monte-Carlo calculations show that switching on the coupling 
between the honeycomb spin planes in a finite size cluster yields a Neel-like AFM spin structure 
with a substantial staggered magnetization at finite temperatures. This may explain the occurrence 
of a robust AFM state in InCu2/3Vi/303 despite an unfavorable effect of structural disorder. 

PACS numbers: 75.50.Ec, 76.30.Fc, 76.60.-k, 75.10.Jm 



In planar honeycomb lattice systems a combination 
of nontrivial topology, strong electronic, spin and or- 
bital correlations and degeneracies yields a rich variety of 
ground states, novel excitations and exotic behaviors that 
currently attract much attention. The recently discov- 
ered exciting phenomena range, e.g., from the quantum 
Hall effect in graphene^, superconductivity in MgB2r 
and intercalated graphite 4 -, to topologically driven quan- 
tum phase transitions in anyonic quantum liquids 5 -. 

Regarding the spin degrees of freedom, an important 
feature of low dimensional spin systems is the presence 
of quantum fluctuations that inhibit long range order of 
the quantum spin-1/2 lattice. Such an effect essentially 
depends on the spin coordination number z. The one 
dimensional (ID) Heisenberg antiferromagnetic (AFM) 
5 = 1/2 chain with z = 2 does not show any magnetic or- 
der even at zero temperature^, whereas long range order 
is possible at T = in the 2D casei as, e.g., in the promi- 
nent S — 1/2 Heisenberg square lattice model with z = 4. 
The honeycomb lattice has the minimum possible coor- 
dination number of any regular 2D lattice z — 3. Thus 
quantum fluctuations there are weaker than in the ID 
case, but stronger than in the 2D square lattice. Hence, 
the AFM order for the honeycomb lattice is fragile^. 

Experimentally, low-D spin systems described by the 
Heisenberg Hamiltonian % = 2 J a f m ^ SiSj are often re- 
alized in structurally three dimensional organic- or tran- 
sition metal oxide (TMO) compounds where strong AFM 
exchange interaction J a f m between unpaired localized 
spins occurs along only one or two spacial directions. Ow- 
ing to residual small 3D exchange couplings such materi- 
als usually exhibit a long range Neel order at a finite tem- 
perature Tn, though, unlike in the 3D magnets, the order- 
ing occurs at much smaller temperatures Tn <C Jafm/^B- 

Recently, a complex TMO compound, I11CU2/3V1/3O3, 



was suggested as a possible candidate for the realization 
of the S = 1/2 honeycomb lattice^ -. In its layered hexag- 
onal structure the Cu 2+ (3d 9 , 5 = 1/2) ions are proposed 
to be arranged in a 2D network of hexagons with the non- 
magnetic V 5+ (3d ) ions in the center of each hexagon. 
The honeycomb layers are separated by sheets of [InOg] 
polyhedra along the hexagonal c axis. The analysis of 
the static susceptibility x(^) nas yielded the exchange 
parameter J a f m ~ 140 K of the honeycomb spin lattice^ -. 
A later structural neutron diffraction study revealed that 
a structural ^{ViCus^} order in the hexagonal planes 
has a finite correlation length £ st ~ 300A and that these 
structural domains are randomly arranged along the c 
axis 1 - 1 -. As has been argued in Ref . [Ill this structural dis- 
order rendering finite values of the spin-spin correlation 
length £ s < £ st concomitant with the substantial spacing 
between the honeycomb planes and with the low spin co- 
ordination number in the planes could preclude a 3D Neel 
order in IJ1CU2/3V1/3O3. Thus the kink-like anomalies 
in %(T)-dependence at 30 K and 38 K previously identi- 
fied with the transition to the long range ordered state 1 - - 
have been tentatively assigned to glass-like order of un- 
saturated spins in domain boundaries 1 - 1 -. 

To obtain insights into the nature of the ground state of 
L1CU2/3V1/3O3 we have investigated its magnetic prop- 
erties using two local spin probe techniques, high field 
electron spin resonance (ESR) and nuclear magnetic res- 
onance (NMR). These experiments were complemented 
by measurements of the field dependence of the static 
magnetization M(B) and by quantum Monte-Carlo cal- 
culations. From ESR and NMR data we have obtained 
evidence for the development of the AFM sublattices in 
the spin system below ~ 20 K which according to the 
M(B) data experience a reorientation transition at a 
field of ~ 5.7 T. Our Monte-Carlo calculations reveal that 
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switching on the interlayer coupling in a finite-size honey- 
comb spin domain yields the formation of the staggered 
sublattices with substantial sublattice magnetization at a 
finite temperature. We conclude that structural disorder 
does not inhibit the Neel ordered state in InCu 2 /3V 1 / 3 03 
at least on a spacial scale of structural domains which, 
owing to their large size, has properties very similar to 
the long range AFM order in the infinite systems. 

Clean stoichiometric powder of InC^/sV^Os was 
synthesized and thoroughly characterized as described in 
Ref. ITD . From it we prepared an 'oriented' sample by 
mixing the powder with cpoxy and letting it harden in 
a magnetic field of 1 T. Owing to the specific anisotropy 
of the Cu 2+ g factor— the c axis of the powder particles 
was oriented perpendicular to the orientation axis (here- 
after o axis) defined by the direction of the applied field, 
as confirmed by x-ray diffraction experiments. ESR was 
measured with a home made spectrometer on the basis 
of the Millimeterwave Vector Network Analyzer from AB 
Millimetre, Paris, combined with a 17 T superconducting 
magnetocryostat from Oxford Instruments Ltd.—. NMR 
data were collected on a Tecmag pulse solid-state NMR 
spectrometer with a 9.2 T superconducting magnet from 
Magnex Scientific. The NMR spectra were obtained by 
measuring the intensity of the Hahn echo versus mag- 
netic field. The T\ relaxation time was measured with 
the method of stimulated echo. Static magnetization was 
measured with a superconducting quantum interference 
device magnetometer from Quantum Design. In all ex- 
periments the magnetic field was applied parallel to the 
o axis, i.e., perpendicular to the hexagonal c axis. 

In the previous ESR experiment performed at a stan- 
dard X-band 'low' frequency v = 9.47 GHz a strong para- 
magnetic Cu 2+ resonance signal clearly visible at high 
temperatures experienced at T < 80 K strong broaden- 
ing and shift and could not be detected below 50 K— . 
Since such a behavior could be indicative of an opening 
of the energy gap for resonance excitations, e.g., due to 
the establishment of the AFM order, we have measured 
in the present work the ESR spectra at low T in a broad 
frequency domain extending into the sub-THz range. 

ESR spectra at various frequencies v collected at T = 
8K for a magnetic field applied parallel to the orien- 
tation axis (B\\o axis) are shown in Fig. HJa). Several 
resonance modes are visible. At high frequencies three 
modes can be detected, Li, L^ and L%. They shift to- 
wards higher fields by increasing the v. At lower fre- 
quencies, the L4 mode can be observed which, in con- 
trast, moves towards lower fields by increasing the v. 
The evolution of L\ - L4 is summarized in the frequency- 
field (y vs. B) diagram in Fig. [TJb). Notably, none of 
the v(B)l\ branches correspond to a paramagnetic res- 
onance condition hi> = gfiBB lcs . Here h is the Planck 
constant and hb is the Bohr magneton. Modes L\, L3 
and L4 are obviously gapped. Extrapolation of L3 and L4 
branches to zero field yields an estimate of the gap value 
A ~ 190 GHz. Mode L4 softens with increasing the field 
strength revealing a critical field B c ~ 5.5 — 6T where 
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FIG. 1: (Color online) Frequency v dependence of the ESR 
signals at T — 8K: (a) Selected spectra at different v. Open 
circles, solid circles, triangles and squares indicate the reso- 
nance field B les of the lines L\ - L4; (b) The v vs. B diagram 
of the resonance modes. The symbols correspond to respec- 
tive Li in panel (a), solid lines are model curves (see the text). 



VL4 0. Above B c mode L2 emerges. In fact, taking the 
above estimate of A and the value of the g factor for B\\o 
axis g = 2.23 obtained from the X-band measurements at 
high T (not shown) it is possible to reasonably model the 
L \ — L4 branches with the AFM resonance (AFMR) rela- 
tions for a uniaxial two-sublattice antiferromagnet^. L3 
and L4 can be represented as hv = A ± g/isi? rcs , where 
the plus and minus signs correspond to L3 and L4, re- 
spectively. Li that appears above B c can be described as 
hv I g\iB = y Bf CS — A 2 . Finally, L\ follows the relation 
hvjg\iB = \/ A 2 + £> 2 CS . Apart from some systematic dis- 
crepancy for the L2 branch the model agrees well with the 
experiment. The behavior of L 2 , L 3 and L4 is typical for 
the situation when a magnetic field is applied parallel to 
some 'easy' axis of an antiferromagnet whereas L\ corre- 
sponds to a 'hard' direction^. The field B c at which for 
L2 and L4 the frequency v — > is the spin-flop field B c 
related to A as B c = (/i/g// fl )A^. Taking A ~ 190 GHz 
one obtains B c ~ 6T. Since B was applied in the honey- 
comb plane of I11CU2/3V1/3O3, there must be therefore 
some 'easy' direction in the plane that defines the AFM 
vector of two sublattices. Obviously the direction of this 
in-plane 'easy' axis of each powder particle is random in 
the sample, because specific in-plane crystallographic di- 
rections in our sample are not defined. Therefore both 
'easy' and 'hard' AFMR modes are present. We note that 
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FIG. 2: (Color online) 114 In NMR line at 35 K and 10 K at 
fields smaller (a) and larger (b) than B c ~ 6T; (c) T^ 1 re- 
laxation rate of 51 V nuclei as a function of B at T = 4.2 K. 
Lines are guides for the eye. Inset: 51 V NMR spectra at fields 
smaller (left) and larger (right) than B c at 55 K and 5 K. 



these modes can be seen at temperatures up to ~ 20 K. 

Further evidence for the occurrence of the AFM sublat- 
tices and related staggered magnetization has been ob- 
tained from the NMR measurements of 51 V (/ = 7/2) 
and 114 In (I = 9/2) nuclei in InCu2/ 3 V 1/3 03. At high 
temperatures the spectra of both nuclei consist of a 
main line and a typical satellite structure arising due 
to the quadrupole interaction. In is positioned outside 
the hexagonal (Cu,V)-plane asymmetrically with respect 
to the Cu honeycomb ring. Hence one could expect at 
low T an uncompensated staggered field at this position. 
Indeed, below 20 K, where AFMR modes are observed 
in the ESR experiment, the 114 In NMR signal experi- 
ences a remarkable change if measured in a magnetic 
field -Bnmr smaller than B c [cf. Figs. 0[a)and (b))]. At 
T = 35 K the line is unsplit regardless the field of the 
measurement. However, at T = 10 K the signal splits in 
two components for Bnmr < B c whereas it remains un- 
split for -Bnmr > B c . This observation strongly supports 
the scenario of compensated collinear AFM sublattices at 
B < B c which turn into the spin flop state at B > B c . 
The 51 V NMR data agree with this scenario. The V site 
is symmetric with respect to the Cu sites. Therefore the 
staggered fields from the honeycomb spin lattice should 
be compensated at this position. In fact, the 51 V NMR 
spectrum does not change qualitatively with decreasing 
the T at -Bnmr < B c and shows the presence of a broad 
background at a high field [Fig.[2jc), inset]. The absence 
of the broadening of the main line suggests that this back- 
ground is due to the influence of structural domain walls 
where uncompensated spins create a distribution of the 
local magnetic fields. Remarkably, though the 51 V NMR 
signal remains unsplit at low T, the 51 V longitudinal nu- 
clear relaxation rate T± shows a substantial field depen- 




FIG. 3: (Color online)Field dependence of the magnetization 
M(B) at different temperatures for B\\ o-axis after a sub- 
traction of a small T-independent spurious ferromagnetic-like 
signal persisting up to room temperature. Inset: (a) - field 
derivative of M(B); (b) T dependence of the reorientation 
field. Error bars correspond to the width of the transition. 



dence: It decreases by a factor of 7 by approaching the B c 
and shows no field dependence at Bnmr > B c [Fig.^c)]. 
An initial rapid drop of Tj - at B ~ 3T could be due 
to a field-suppression of magnetic fluctuations caused by 
impurity free-like spins. However, according to the mag- 
netization measurements (see Fig. [3] and its discussion 
below) these spins are fully polarized at stronger fields. 
Thus a further substantial reduction of T± is likely to 
be related with the approach to the critical field B c . The 
sensitivity of T^ 1 to the spin structure may be related to 
magnetic anisotropy that defines a specific 'easy' orienta- 
tion in the honeycomb spin plane and possibly produces 
a small uncompensated fluctuating field at the V nuclei 
in the collinear phase that become ineffective for the T\ 
relaxation after the spin reorientation has occurred. 

Motivated by the findings of the local ESR and 
NMR techniques we have performed the magnetiza- 
tion measurements on an oriented powder sample of 
InCu 2 /3Vi/303 for a magnetic field applied parallel to 
the o axis. The temperature dependence of the magnetic 
susceptibility x = M/B at B = 1 T demonstrates the 
same behavior as observed before in Ref. [ll| (not shown). 
The field dependence of the magnetization M{B) at dif- 
ferent T is presented in Fig. [3J The curvature of the 
M(B) at low fields is likely due to a small amount 
of free-like spins, which also yield a Curie-like tail in 
the T dependence of \ (Ref. [lj). The M(B) depen- 
dence at the lowest temperature T = 2 K exhibits a pro- 
nounced step-like increase of M at B ~ 5.7 T. It can be 
straightforwardly interpreted as a critical field B c for the 
spin reorientation, in a nice agreement with ESR and 
NMR data. By increasing the T the step continuously 
broadens, slightly shifts towards higher fields and finally 
smoothly vanishes above ~ 20 K [Fig. 03 insets (a,b)]. 

Thus the results of three different techniques strongly 
suggest the formation of collinear AFM Neel sublattices 
in the S = 1/2 honeycomb plane in InCu2/3Vi/303 at a 
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FIG. 4: T dependence of the staggered magnetization m a t 
of the honeycomb lattice with interplane couplings J z = 
0.1, 0.2, 0.3 and 0.4J a f m for the system size 14 x 28 x 14. Inset: 
m st vs. T/J a f m for J z — 0.2J a f m for different system sizes. 

finite temperature. As discussed in Ref. [Ill the absence 
of the magnetic anomaly in the specific heat C p does not 
exclude the occurrence of the magnetic phase transition 
in a low-dimensional system such as InCu 2 /3Vi/3 03 due 
to residual interactions in the third dimension. If the 
ordering temperature Tn is substantially smaller than 
Jnfm/kB, most of the magnetic entropy S'mag is already 
gone at T ~ Jafm/^B due to the 2D in-plane AFM corre- 
lations. Thus an additional change of S'mag at Tn might 
be small and not visible in the C p owing to a substan- 
tial phononic background. The question arises if this 3D 
situation is realized in InCu 2 /3Vi/3 03 despite the occur- 
rence of disordered finite size structural domains. To 
address this issue we have performed a computational 
study of a finite size spin- 1/2 clusters on the honey- 
comb lattice by means of a continuous Euclidean time 
Quantum-Monte-Carlo algorithm that includes the cou- 
pling in the third dimension-^. We find that the Neel 
type structure and the staggered magnetization m st , oth- 
erwise present in the 2D case only at T = 0, develops in 
a rather narrow interval at a finite temperature that de- 
pends on the strength of the interlayer coupling J z . As 
shown in the main panel of Fig. [4] a temperature inter- 
val where m s t(T) exhibits a step-like increase is mainly 
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determined by the value of J z whereas the sharpness of 
this step increases with the increase of the cluster size 
(inset of Fig. QJ. Note that in InCi^/^Vi/^Os the size 
of in-plane structural domains amounts to ~ 300A (Ref. 
[TTh which corresponds to ~ 90 x 90 spin sites, i.e. a much 
larger number than in our model calculation. Thus one 
should expect quite a sharp set up of the two-sublatticc 
collinear AFM Neel spin structure and the occurrence 
of the staggered magnetization in I11CU2/3V1/3O3 that 
experimentally would be difficult to discriminate from 
a true magnetic phase transition expected in the ther- 
modynamic limit. Finally, regarding modelling of the 
finite size systems, we mention here that a calculation in 
the framework of the (anisotropic) Heisenberg model of 
the ESR response of the AFM cluster comprising 8 spins 
only already yields AFMR modes of a substantial inten- 
sity that for the 'easy' axis direction have qualitatively 
similar properties as those shown in Fig. QJ 5 . 

In summary, by studying ESR, NMR and static mag- 
netization of the oriented powder sample of the honey- 
comb spin-1/2 lattice compound I11CU2/3V1/3O3 we have 
found strong experimental evidence for the formation of 
the Neel type collinear AFM spin structure at tempera- 
tures below ~ 20 K and a respective development of the 
staggered magnetization. A reorientation of spin sublat- 
tices in a field of ~ 5.7 T has been clearly identified in the 
magnetization and magnetic resonance data. The AFM 
ordering of spins in InC^/sV^Cs takes place despite 
the low spin coordination number, the pronounced struc- 
tural two dimensionality as well as the presence of the fi- 
nite size incoherent structural domains. Our experimen- 
tal results are corroborated by the Quantum Monte-Carlo 
study of the finite size honeycomb spin layers coupled in 
the third dimension. It reveals a sharp development of 
the staggered magnetization at a finite temperature that 
in the thermodynamic limit would correspond to a phase 
transition to the AFM long-range ordered state. 
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